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There is little consensus on the nature of the glass state and its relationship to other strain states in ferroelastic 
materials which show the shape memory effect and superelasticity. We provide a thermodynamic interpretation 
of the known strain states, including precursory tweed and strain glass, by mapping the problem onto a spin 
model and analytically obtaining the phase diagram using real- space renormalization group methods. We further 
predict a spontaneous transition from the glass state to the ordered martensite phase. We verify this prediction by 
mapping out the experimental phase diagram for the ternary ferroelastic alloy Ti5o(Pd5o-xCrx) and demonstrate 
the emergence of the spontaneous transition. Our work thus provides a consistent framework to understand the 
various experimental and theoretical studies on the glassy behavior associated with ferroelastic materials. 
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I. INTRODUCTION 

Ferroelastic materials undergo first order transitions that 
are characterized by a lattice strain or shuffle, and trans- 
form from a high temperature strain-disordered paraelastic 
state (austenite) to a low temperature strain-ordered ferroe- 
lastic state (martensite) where long-range elastic interactions 
are important. It is known that the thermodynamics of such 
phase transitions is strongly influenced by the presence of 
disorde^JEl^ In particular, statistical compositional fluctua- 
tions play a fundamental role in bringing about a precursory 
strain state known as tweed. This is a crosshatched pattern 
observed well above the martensitic-transformation start tem- 
perature^^'^. Recent experiments on ferroelastic alloys have 
shown that by introducing disorder via doping point defects or 
compositional variations beyond a critical value, an abnormal 
glass-like state which is a frozen state of local strain order, 
can be generated below a transition temperature^^ This 
glass phase was initially observed in Ni-rich Tiso-a^Niso+x 
where the austenitic B2 parent structure appeared to persist to 
K above a compositional threshold x > 1.3, below which 
it transformed to a martensitic B19' phase^^ This so called 
"strain glass" is of interest not only from a theoretical point 
of view but it has been shown to display superelasticity and 
shape memory effects, which are typically seen in austenite 
and martensitic states 

Experiments and theory have so far provided little under- 
standing of the nature of the glassy behavior in ferroelas- 
tic materials and its relationship to other strain states. For 
example, the tweed has been the subject of a number of 
theoretical studies and has variously been interpreted as a 
"glass" phase^ Although recent experiments distinguish 
strain glass from tweed^^, the claims are largely based on di- 
agnostics that monitor certain static and dynamic aspects of 



the materials; these include the broken ergodicity of static 
properties (e.g., the strain measured by zero field and field 
cooling (FC-ZFC)), or frequency dispersion of dynamic prop- 
erties^^. However, non glassy phases, such as polytwins in 
martensite or nano structured ferromagnetics (e.g. FeReCr), 
also show similar non-ergodic static behavior and frequency 
dispersion In addition, although numerical simulations 
based on continuum Landau descriptions in the presence of 
disorders can reproduce some form of the experimental re- 
sults such solutions are also based on similar empirical 
diagnostics and tend not to be predictive. Therefore, there 
is a need for a predictive approach using analytical techniques 
that would allow the various ordered and disordered ferroelas- 
tic states to be distinguished and which can be experimentally 
verified. 

In this paper, we provide a thermodynamic interpretation 
of the known strain states in a ferroelastic system, including 
austenite, martensite, tweed and strain glass. After mapping 
the problem onto a spin model, we use real space renormal- 
ization group (RG) methods to analytically calculate thermo- 
dynamic phases in terms of RG attractive fixed points. The 
RG approach progressively integrates out microscopic de- 
grees of freedom so that the attractive basins characterize the 
physics at large scales. The values of the interaction strength 
and strength of disorder uniquely characterize the different 
phases, including the frustrated glassy state, thereby allow- 
ing the phase diagram in terms of temperature and disorder 
to be determined. We thus theoretically determine the phase 
diagram of a model ferroelastic, and the phase diagram pre- 
dicts a spontaneous transition from strain glass to marten- 
site if the strength of disorder in the system has intermedi- 
ate values. This latter aspect was not recognized previously, 
although the spontaneous transition has been recently spec- 
ulated^^. We verify these predictions by mapping out the ex- 
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perimental diagram for a ternary ferroelastic Ti5o(Pd5o-a;Cra;) 
alloy (where x is the atomic concentration of Cr) and exam- 
ining its consequences. Moreover, our calculations show that 
the tweed phase can be interpreted as a thermodynamic equi- 
librium phase. Our model thus provides a unified, consistent 
framework to understand the various experimental and theo- 
retical studies on the glassy behavior associated with ferroe- 
lastic materials. Moreover, we find that the long range elastic 
interaction or the precise form of the disorder are not crucial 
to describing the phase diagram. These ideas are well recog- 
nized for spin glasses but are not as well known within the 
strain glass community. 

The plan of our paper is as follows. In the first section, 
we review how effective spin models can be deduced from 
Ginzburg-Landau (GL) functionals, and we discuss the vari- 
ous means to include quenched disorder in the model. The 
second section is devoted to the calculation of the analyti- 
cal phase diagram from our model in the presence of dis- 
order. We explain how a real-space RG procedure can be 
implemented and we present several approaches and approx- 
imations to solve the resulting equations. The results are 
compared to Monte-Carlo simulations and we provide a dis- 
cussion of the relevance of simplified numerical algorithms 
- be it a one- spin flip Monte-Carlo or the steepest-descent 
method for continuum GL models- for these types of prob- 
lems. The role of the long-range interactions is also discussed 
in relation to previous work. Section |IV| contains our exper- 
imental results on TisoCPdso-^Cra^) alloys. We measure the 
transformation behavior of the alloys by means of dynamical 
mechanical analysis (DMA) and the electric resistivity. In- 
situ synchrotron X-ray and transmission electron microscopy 
(TEM) were employed to detect the predicted spontaneous 
phase transformation. Our experiments lead to a phase dia- 
gram that is very similar to the one we establish analytically. 
Finally, we discuss our results in the broader context of ferroe- 
lastic transitions in general and the possible generalizations of 
our theoretical work to more complex transitions. 



II. PSEUDO-SPIN MODEL AND QUENCHED 
IMPURITIES 



Our approach is to use a model that captures the salient 
physics and micro structure associated with ferroelastic trans- 
formations. We will use a pseudo-spin model derived from a 
continuum formulat ion b ased on a Landau potential and strain 
compatibility forces The addition of quenched disorder 
then permits studies of the model using the tools of statis- 
tical mechanics. We review in this section the crucial steps 
in the derivation of such spin models starting from standard 
Ginzburg-Landau (GL) free energy functionals. We then dis- 
cuss how the effect of impurities in the original GL theories 
can be taken into account in spin models. 



A. Spin model 

The idea of using a spin model for ferroelastic transitions 
was introduced in Ref. This idea was later generalized 
to many other transitions and shown to capture the salient 
physics of continuum GL modelJ^. For simplicity, we will 
study here 2D square to rectangle (SR) transformation driven 
by the deviatoric strain 62 = ^(en — 622), where e^^^ are 

components of the strain tensor e^jy = ^ ( + |^ ) defined 
in terms of displacements, u. The generalization to more com- 
plicated - and more realistic - transitions will be discussed in 
Sec, [vl The Landau free energy of the system is Fql [^2] = 
Eq jdPr [/l + /g + Ilr], where the local free energy for 
the first-order transition is = (r — V)e^ ^ e\(e\ — \Y , 
the gradient term fc — ^ |Ve2|^, accounts for the cost 
of creating interfaces between different variants and /lr = 
^ J d?r' e2{r)U{f— 7^)62(7^) is the elastic long-range force. 
The scaled temperature r = is expressed in terms of 

the transition temperature Tg^ and is the temperature of the 
austenite stability limit. The long-range term comes from the 

so-called compatibility equation V x x = on the 

strain tensor e, that ensures that the displacement field is sin- 
gle valued. This leads to the anisotropic Kernel U{f — f), 
which reads in Fourier space 



U{k) 



(1) 



where ^(r) = [| (l + ^1 - ^) 



The factor v{k) = 1 — Sk^Q ensures that this long-range term 
vanishes for a uniform strain field. In all what follows, the 
ratio = ^ will be kept constant. Back into real space 

we roughly have U 

The strain 62 is the order parameter for this transition. The 
Landau term = (r — 1)62+62(62 — 1)^ ensures that the sys- 
tem undergoes a first-order phase transition at r = 1. When 
r < 4/3, this term has three minima, which correspond to 

phases called austenite (62 = 0) and martensite (62 = ±£:(r)), 

1/2 

. For 1 < r < 4/3, 
the austenite is a metastable state. To obtain non-uniform 
textures, one minimizes the GL functional Fql [^2] as usual. 
Note that at this stage the GL continuum model is in itself 
a mean-field description of the microscopic problem, and as 
such it does not contain any information about fluctuations 
around the mean-field solution given by the minimization of 
^GL [62]. An alternative point of view is to consider a for- 
mal partition function 2^ = J V[e2]ex.p (— /3F[62]), where the 
strain configurations are summed over with a weight given by 
the usual Boltzmann weight. The Laudau (or Mean-field) ap- 
proximation is then nothing but a saddle point calculation of 
this functional integral. One can simplify the calculation of 
Z by retaining only the strain minima given by the Landau 
part of F, transforming then 62 into a discrete variable that 
we rewrite as 62 (r) = £{r)S{f) with the pseudospin S{f) = 
0, ±1. This is the main point of the pseudo-spin approxima- 
tion, instead of keeping only one configuration which mini- 
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mizes the full GL functional Fql [^2], we retain in the par- 
tition function all the configurations that minimize the local 
Landau term /l = (r— 1)62+62(62 — 1)^. Within this approx- 
imation, the partition function reads Z XI {5} ^^P (~/^^)' 
where the Hamiltonian H is obtained from the continuum the- 
ory as H = F ql [e2{r) £{r)S{f)]. After a proper dis- 
cretizatiorP^I^, it reads 



= -J{t\ 



E 

<i,3> 



5,5, + A(r) ^ 5^ + ^ ^ S,U,,S„ 

i ij 

(2) 

where A(t) = Do(t)/2 (^^(t) + 4^^) and J(t) = 
Do(T)e2,withDo = 2l3Eoe{Tf,gL = (T-l) + (£2(r)-l)2. 
The states Si = ±1 correspond to the two rectangular variants 
while Si = represents to the square austenite. Note that this 
spin model only makes sense for r < 4/3, as for r > 4/3, 
only the state S = is allowed. We have therefore mapped 
our continuum GL theory onto a lattice spin— 1 model with 
Hamiltonian ([2]). This Hamiltonian with Ai = is known in 
the spin literature as the Blume-Capel model^^ and it repro- 
duces all the well-known features of the transitiorP^, in par- 
ticular, it has a first-order transition around r o:^ 1. 



There are several points that are worth mentioning at this 
stage. It is usual to deduce GL theories from lattice mod- 
els, not the other way around. The reason for this is that GL 
theories are easier to deal with and contain interesting coarse- 
grained, mean- field information about the initial microscopic 
model. In our case, we chose to start from a simplified GL 
theory to map onto a spin-1 statistical model that may seem 
harder to handle. On the other hand, spin models will turn 
out to be very convenient when introducing quenched dis- 
order as connections to usual spin glass models can then be 
mad^El; It is also important to realize that when progressing 
from mean-field GL theory to a classical spin model we have 
actually put in additional information about thermal fluctua- 
tions, as the former a priori does not contain any information 
about fluctuations around the mean-field solution. The Hamil- 
tonian ^ is an effective model, with temperature-dependent 
coefficients and is not an accurate microscopic description of 
the phenomenon that we are studying. The additional infor- 
mation contained in the effective Hamiltonian ^ is somehow 
arbitrary and other choices of reasonable Hamiltonian would 
be possible. For example, one could as well have introduced 
two different temperatures, the physical temperature T that 
would appear in r, and another more artificial temperature Tg 
that would appear in the Boltzmann factor exp(— i^/Tg). As 
was suggested in Ref.^, one can then also let Tg = so that 
one would need to minimize H in order to obtain statistical 
properties. We choose to introduce somewhat artificially ther- 
mal fluctuations by taking T = Tg, a point of view completely 
analogous to what is usually done with the well known 0^ the- 
ory. This is a crucial point when dealing with such effective 
spin models. 



B. Spin Model with Quenched Disorder 

Now that we have seen how to model ferroelastics using 
spin models, we consider the influence of impurities in the 
model. Our aim is to understand how quenched impurities 
affect the phase diagram of ferroelastic materials. This is 
usually done phenomenologically in Laudau theor>^^ by in- 
troducing a random field or by introducing randomness in 
the transition temperature Tg with non-zero spatial correla- 
tions^^ In this paper, we follow the idea that one need 
not worry about the precise microscopic form of the disor- 
der induced by impurities as one expects the effect of quench 
disorder in the interactions of the model to yield somewhat 
universal features. Recall that our spin model is an effec- 
tive model, so the results we are after are generic features of 
the phase diagram, such as its topology, rather than an ac- 
curate, quantitative description to be compared directly with 
experiments. Note that this train of thought is very similar to 
what was done in spin glass theory for usual ferromagnets (see 
^.g.EEIl). The point is then to start from the pure (disorder- 
free) Hamiltonian ^ and to take the nearest-neighbor cou- 
plings to be quench independent random variables Jij , drawn 
from the distribution V{Jij), with mean J(r) and with vari- 
ance a J. As we will see later, the precise form of this distribu- 
tion is irrelevant to the global topology of the phase diagram. 
This is a strong argument in favor of the "universality" dis- 
cussed previously. We believe that this way of introducing dis- 
order tends to be more satisfying than the very specific, some- 
times fine-tuned, methods usually used in the literature. The 
parameter aj can be thought of as a measure of the quenched 
disorder in the system, in particular, for a j = we recover a 
pure system. For future reference, we give the explicit form 
of our disordered Hamiltonian 

I3H = -J2 J^Ar)S^S,+A{r) ^ Sf+^ ^ S,U,,Sj. 

<i,j> i ij 

(3) 

where the couplings Jij are quenched variables drawn from 
the distribution V{Jij). In this paper we will use two different 
types of distributions 



nJzj) = l^iJzj - Ji) + \^{Jz3 - (4a) 



njzj) 



1 ^^^( iJ,-Jir)y 



(4b) 



/27raj ' V 

with Ji = J[r) -\- aj and J2 = J(r) — crj. In the following 



the distribution (4a) will be referred to as "Bimodal" whereas 
we will denote (4b) as Gaussian. These equations define the 
model we will attempt to solve in the next section. 



III. PHASE DIAGRAM 

We now consider how real-space Renormalization-Group 
(RG) can be used to compute the full phase diagram of the 
the spin model ([5]). We also present preliminary Monte Carlo 
results and discuss their relationship to the phase diagram of 
the model. 
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A. Real-Space Renormalization Group for the pure model to the scheme 



Our goal in this section is to obtain an analytic phase di- 
agram of the model defined by eq. ^ in the limit Ai = 0. 
From now on we will only take the limit of vanishing long- 
range interactions Ai = as we do not believe they are im- 
portant to understand the phase diagram of ferroelastic mate- 
rials. We will return to this later. There are several ways to 
deal with a Hamiltonian such as ([2]), an obvious straightfor- 
ward method would be to perform a mean- field approxima- 
tion. This yields (see e.g. Ref.^^ for a description in the con- 
text of martensitic transitions) a first order phase transition at 
T 1 between a martensite phase with m = (S) ^ and 
an high-temperature austenite phase characterized by m = 0. 
Here we choose a different, usually more reliable method: the 
real- space Renormalization Group (RG). Its application to our 
pure spin-1 model goes back to Ref.^^, and is of course well 
known. We review here some of the essential elements before 
proceeding to the case of disorder. 

The idea is to approximate the square lattice by a hierarchi- 
cal lattice, obtained by iterating the construction scheme 




(5) 

Several other choices of hierarchical lattices to approximate 
the square lattice are possible, however the one we choose is 
known to be a very good approximation in the case of spin- 
1 models (see e.g. Ref.^^). We define our spin model on 
such a lattice by the Hamiltonian I3H = —J^^i SiSj + 

^Yj<ij>i'^i~^'^j)' Note that we expressed the crystal field 
term A as a sum over the links instead of a sum over sites 
in order to take into account the fact that the coordination 
number of the sites in the lattice ^ is not constant. It is 
well-known that such precautions are needed when dealing 
with the renormalization of local field terms such sls —h^Si 
or specifically here ^^Sf. This corresponds to assigning 
weights to the sites in the crystal-field interaction term (A) 
according to their coordination numbers which are, unlike the 
square lattice, not uniform. This inhomogeneity is necessary 
so that the square lattice is correctly approximated using a hi- 
erarchical lattice. 

On such a lattice, the partition function can be computed 
exactly by summing progressively over the spins. The renor- 
malization procedure consists in decimating spins according 




(6) 

so that the renormalization process corresponds to inverting 
the arrows in eq.([5j. More formally, if one denotes by i^i234 
the Hamiltonian of the left-hand side of eq ([6]) and by its 
right-hand side, one has 

exp [-pH[2] = Tr exp [-pH.^s^] , (7) 

where Tr amounts for a partial summation performed on spins 
5'3 and ^4. This equation then gives the relationship between 
i^i234 and the renormalized Hamiltonian H[2,^o that one can 
draw the evolution of the coupling parameters by repeating 
this decimation procedure. It is not hard to see that this trans- 
formation is exact if and only if one introduces an additional 
coupling —K'^^.j^SfSj that is generated upon renormal- 
ization. Although this biquadratic interaction is absent in 
our original model, it is generated by the RG procedure, and 
must be taken into account to follow the exact RG flow. In 
our case, the elementary Hamiltonian H1234 reads = 
- E [JSiSj^KS,^Sj^]^lA E [5,2 + 5/], where the 

sum are performed over the links ((13), (14), (42), (32), (34)) 
in eq. ([6]). It is worth emphasizing again that the additional 
coupling term —K^^.js^Si'^Sj'^ must be introduced in the 
Hamiltonian, otherwise there is no analytical solution to be 
found for eq. ([7]). The renormalized Hamiltonian then 

reads /3i7{2 = -J' 8182- Si^ 82^ ^ ^ (^i' + ^2')+C^ 
where C becomes a simple multiplying constant exp(— 
that contributes to the renormalization of the free energy, ir- 
relevant for our purposes. 

It is straightforward to solve eq. ^ for all 81,82 • This 
yields a non-linear map {J'^A'^K') = 7^ [ J, A, i^T] that gives 
the exact renormalized couplings of the effective Hamiltonian 
in terms of the initial ones. Iterating this map, one can 
obtain the effective Hamiltonian describing the physics of the 
system at large scale, where the microscopic degrees of free- 
dom of the system have been summed over. Basins of at- 
traction of this map correspond to thermodynamical phases, 
characterized by attractive fixed points. For example, the 
martensite phase (ferromagnetic phase in the spin language) 
has J* = +00 and A* = —00, so at large scales the ef- 
fective Hamiltonian describing the system forces the spin to 
be in the 5* = ±1 states, with a strong correlations between 
nearest neighbors. Meanwhile, austenite can be described by 
the attractive fixed point J* = and A* = 00 so it cor- 
responds to a disordered phase that favors 5* = 0. To ob- 
tain the phase diagram of our spin model ([2]) for A 1 = 0, 
we simply iterate the map 1Z starting from the initial point 
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(J, A, K) = (J(r), A(r), 0). We find that the system flows 
to the martensite fixed point for r < 0.96, and to the austenite 
fixed point for r > 0.96. Moreover, the study of the largest 
eigenvalue of the linearized Renormalization- Group matrix 
around the transition at r 0.9 6 can be used to prove that 
the transition is of first ordei^SSI^ as expected. 



Phase 


OP characterization 


RG fixed point 


Austenite 
Martensite 

Tweed 
Strain glass 


m — q — small 
m/0,p/ 0,^7^0 

m — q — large 
m = 0,p/ 0,^/0 


A* = +00, J* = 0, cr} = 

A* = — oo, J* = oo, ^ = 
A* = -oo, J* = 0, cr} = 
A* = -oo, J* = 0, cr} = oo 



B. Real-space RG with Disorder 

The idea of using spin glass models to describe disor- 
dered ferroelastic materials originates from Refs.^^^. The 
replica/Mean-Field approach of this model^^ yields a strain 
glass phase but also several unwanted features that are con- 
trary to both experiments and Monte-Carlo simulations. The 
absence of a tweed precursor is an example. We use here a 
different approach relying on real- space RG which we believe 
is simpler and usually more reliable. Although our starting 
point is the same Hamiltonian ([3]), we will see in this paper 
that RG provides interesting new results that can be tested ex- 
perimentally. 

The implementation of the RG procedure for disordered 



systems is very similar to the pure case in Sec. Ill A The 
real- space RG in disordered systems still relies on a decima- 
tion procedure to progressively sum over the spin variables, 
integratin g out the microscopic degrees of freedom (see 
e.g. RefsP^ISEIl). After each summation step, the effective 
Hamiltonian describing the remaining degrees of freedom 
is assumed to be the same as the original one but with 
renormalized coupling coefficients. This turns out to be exact 
for the hierarchical lattice ([5]). The only difference is that for 
disordered systems the RG allows us to follow the probability 
distribution of the couplings upon the renormalization proce- 
dure instead of the couplings themselves. In order to do so, 
the main point is to obtain the RG equations in the case of a 
completely inhomogeneous system. Even though the initial 
Hamiltonian (jSj (with Ai = 0) has pure K = and A inter- 
actions, the RG flow will lead to randomness in all couplings, 
so one has to consider the more general Hamiltonian PH = 



TABLE L Characterization of the different thermodynamic phases 
in terms of the calculated Renormalization Group (RG) fixed points. 
The austenite and tweed phases are clearly separated because one 
fixed point has value A* — +oo (so the state = is favored) 
whereas the other has A* = — oo. We also show the corresponding 

OP values m — (Si), p = (Sf), and q = {Si)"^ for the phases as 
characterized in mean field and replica theory. 



Kij' = 7^[{Kij}] then gives the 4 renormalized couplings 

^ij' — (^125^125^127^12) ^ function of the 20 initial 
ones {Kij} = (K13, K14, K34, K23, K42), where we used 
the labeling of eq. ([6]). The RG recursion relation for the joint 
distribution 7^ (Kij) = V{Jij, K^j, Aij, A]-) then reads 



E 

(id) 



The initial condition for the RG flow is the Hamiltonian ^ 
with Ai = 0, so that before renormalization, K^j — 0, 
Aij = A(r) and A^ = 0. Remark that we once again 
considered the crystal field term as living on the edges of the 
lattice, and not on the sites. This should be reminiscent of 



our treatment of the pure case in Sec. Ill A The solution of 
eq. ^ in this case yields a map IZ that gives the renormalized 
couplings (7(2, ^12 7 ^12 5 ^12) terms of the original ones 

( Jl4, Ji3, . . . , Ki4^ i^i3, . . . , Ai4, Ai3, . . . , AJ4, AJ3, . . . ). 

One then analyzes the evolution of the joint probability distri- 
bution V{Jij, Kij^ Aij^ aIj) upon renormalization, starting 
from the initial distributions V{Aij) = S{Aij — A(r)), 
P(A|^.) = 5{Ay, V{K,j) = 5{Kij) and P{Jij) given 
by either ( |4a| or ( |4b| ). Let us now be slightly more 
explicit. Let Kij = ( J^j, i^T^j, A^j, A^) denote the 
four couplings living on the link {ij). The RG maps 



P'(Ky') 



ndKyP(Kij) 



5 (Kij' - 7^[{Ky}]) 



(8) 

where the product is over the five links of eq. ([6]), so there 
are 20 integrations to be performed at each RG step. We em- 
phasize again that even though one starts from pure K and 
A couplings, randomness in these parameters will be gener- 
ated by the RG procedure. Similarly, even if one starts from 
a bimodal Gaussian distribution for V{Jij), the resulting dis- 
tribution obtained after several renormalization steps could in 
principle be much more complicated. 



C. Basins of attraction and thermodynamical phases 

Before solving explicitly eq. ([5]), let us discuss how the 
phases are characterized within this formalism. The RG 
distribution V{Jij^Kij^ Aij^ A\-) will typically flow to at- 
tractive fixed points that will characterize thermodynami- 
cal phases. In principle, one would need an infinite num- 
ber of parameters to characterize RG-invariant distributions 
V* {Jij , Kij , Aij , A^ ) . However, for our purposes it turns 
out that the thermodynamic phases can obtained as RG basins 
of attraction characterized only by the values of A* and 
(J*,cr}) at the fixed point. Here A* and J* are the mean 
values of J and A at the fixed point, while a} is the standard 
deviation of J. 

The martensite (ferroelastic) phase corresponds to the usual 
ordered ferromagnetic phase (J* = 00, A* = —00, a}/ J* = 
0). In terms of order parameters (OPs) used in mean field 
theory and replica calculations (see e.g!^, this phase is char- 
acterized by a non-zero magnetization m = {Si) 7^ 0, where 
the bar represents average over the disorder and the brack- 
ets average with respect to Boltzmann weights. We find two 
paraelastic disordered phases with J* = 0, a} = 0, and 
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A* = ±00. The case A* = 00 corresponds to the austen- 
ite case as it favors S = whereas the case A* = —00 is 
interpreted as a disordered phase of martensite clusters which 
we identify as the tweed phase. The OP that allows to dis- 
tinguish between both phases is the "martensite volume frac- 
tion" p = (Sf). We also remark that the tweed precursor 
we find is ergodic and non-glassy, consistent with recent ex- 
periments^^. This is to be compared with the hypothesis of 



0^ that interpreted tweed as a glassy phase. Our model 
results are that tweed is a thermodynamic phase, rather than a 
metastable precursor. Note also that this tweed phase is not 
captured by a mean-field/replica-symmetric analysis of our 
2D model^^. The last phase that we encounter corresponds 
to a spin/strain glass with infinite randomness (a} = 00) and 
(J*=0,A*=oo). The effective Hamiltonian describing the 
system at large scales has features in which variants S = ±1 
are favored (because A* = 00), and the values J* = 
and cr} = 00 imply random infinite couplings on each bond, 
thus denoting frustration. This phase is also characterized by 

the Edwards-Anderson order parameter q = (Si)'^ which in 
the replica language corresponds to the overlap between two 
replicas q = (^SjSf) of the systenP'. The characterization in 
terms of RG fixed points or OPs of these four thermodynamic 
phases is gathered in Tab. |T| 
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FIG. 1. Analytical phase diagrams from the exact numerical resolu- 
tion of the RG equations on the hierarchical lattice. Main figure: the 
initial distribution of the disorder was taken to be bimodal . Inset: 
initial Gaussian distribution. Notice the absence of spin/strain glass 
phase (see text). 



D. Analytical phase diagram 

1. Numerical Resolution 

We now discuss how one can solve ([8]). An obvious 
method would be to sample numerically the joint distribu- 
tion V{Jij , Kij , Aij , A J ) . This allows one to solve explicitly 



eq ^ numerically. We used pools of 160000 values to sam- 
ple the distribution V{Jij, K^j, Aij^ aIj). Results are shown 
in Fig.[T]for both bimodal ( |4a| and Gaussian (|4b]) initial distri- 
butions. Parameters appropriate for our example are £^0 = 3, 
^ = 0.5, Teq = 1, and Tc = 0.9. Note that strictly speaking, 
our spin model is not defined for r > 4/3, as the spin approx- 
imation only yields one state 6* = in that case. However, we 
can still think of this region as being in the austenite phase, 
as this is the only phase allowed by the spin approximation. 
One of the most striking features of these phase diagrams is 
the appearance of an intermediate tweed phase between the 
austenite and martensite phases as the disorder is turned on. 
This is consistent with what is known experimentally (see ex- 
perimental section). Meanwhile, the phase diagram has the 
same topology for both types of disorder. This is a strong re- 
sult in favor of the "universality" discussed in Sec. ( |II A| ). 

One notes the absence of the spin/strain glass phase, which 
is to be expected. Strictly speaking, the model ([3j in two di- 
mensions should not have a spin-glass phase at finite temper- 
ature^^ . To be more precise, the Hamiltonian ([Sj is believed to 
have a spin glass lower-critical dimension dc lying somewhere 
between 2 and 3. This means that for d < dc (low dimension), 
and in particular for d = 2, the spin glass phase should be de- 
stroyed by thermal fluctuations. Meanwhile, it was recently 
proved^ ^ that the very same model (forgetting about the tem- 
perature dependence of the coefficients) has a spin glass phase 
at finite temperature in 3D. (Note that the SR transition em- 
bedded in 3D - or our spin model in 3D - corresponds to a 
slightly constrained tetragonal to orthorhombic transition). Of 
course, a spin glass phase can also be found within a mean- 
field/replica solution of the modePI, as mean-field effectively 
corresponds io d = 00. We expect this spin glass phase to 
reappear on higher-dimensional hierarchical lattices^. Al- 
though one might think that the absence of spin glass phase 
in 2D could be of crucial importance experimentally, this is 
actually irrelevant to us, for the following reasons: 

• Even though a spin glass phase should not exist in 2D, 
we still expect kinetic features to make the system look 
"glassy". This will be discussed in more detail in the 
Monte Carlo section below. The point is that the dis- 
tinction in experiments and numerics between a gen- 
uine spin glass phase and a glassy kinematic behavior 
is actually a very subtle issue. 

• As our spin model was derived from a mean-field Lan- 
dau energy without thermal fluctuations, we are inter- 
ested in generic (mean-field like) features of our model. 
This is a consequence of the discussion at the end of 
Sec. ( [II A) . The spin glass lower critical dimension is 
not a meaningful quantity in our case as our model was 
derived from a mean-field Landau energy in the first 
place. Therefore, the fact that the somewhat artificial 
fluctuations - introduced in going from the continuum 
GL theory to our spin model- may or may not destroy 
the spin glass phase is not relevant , as these fluctuations 
do not exist in the original Landau model. Our model is 
thus only meaningful in a mean- field-like context. 

Therefore, one need not worry about this issue of lower crit- 



7 



ical dimension here, it would be meaningful only if our spin 
model were a precise microscopic description of ferroelastics. 



2. Projection Approximation 

Although the direct numerical resolution of eq. ([8]) is the 
most straightforward way to proceed, it is also useful to have 
an approximate way of solving this system, thereby allowing 
one to perform the calculations analytically. A possibility to 
avoid this rather cumbersome numerical procedure is to make 
a further approximation in the case of Bimodal disorder (|4a|). 
One can remark that rather than follow the full evolution of 
these distributions, one could enforce the renormalized dis- 
tributions to be the same as the initial ones but with renor- 
malized parameters. That is, we enforce the distribution of 
the Jij couplings to remain Bimodal upon renormalization, 
keeping K and A constant. Although this approximation may 
appear too drastic, it is quite common in spin glass related 
problems and we will show in the following that it captures all 
the important features of the phase diagram of the model ([3]). 
Note that such an approximation typically yields results char- 
acteristic of higher-dimensional systems (see e.g. Ref.l^and 
references therein), and thus should alter the spin glass lower 
critical dimension. As discussed previously, this is not impor- 
tant for our purposes. In the following, we will refer to this 
approximation as the "projection approximation", as it indeed 
consists in projecting the renormalized distributions onto the 
initial ones. 

Fig.|2ja) shows the phase diagram in the plane (r, aj) ob- 
tained by iterating the RG map in the projection approxima- 
tion case. In the absence of disorder (aj = 0), we find a 
first-order phase transition between the austenite and marten- 
site phase with r 1, as expected. As one increases the dis- 
order, an intermediate tweed phase arises before it transforms 
into low temperature phase (either martensite or strain glass). 
For large enough disorder and low temperature, we find a spin 
glass phase that we interpret as a strain glass. Interestingly, 
when the disorder of the system is in the intermediate regime 
(1.3 < (jj < 2.3 in our model), we find there should exist a 
spontaneous phase transition from strain glass to martensite, 
this is a predictio n that wa s not obtained by previous numeri- 
cal simulation^ 13|19|23|24|33| based on Landau theory. The exis- 
tence of a spin glass phase in this calculation is related to our 
projection approximation, which is known to produce results 
characteristic of higher-dimensional systems (it is therefore 
legitimate to think of this approximation as a kind of mean 
field). Note also that except for the spin-glass phase, the exact 
numerical resolutions of Fig. [T] and the projection approxima- 
tion results of Fig.[2ja) are very much alike. 

E. Monte Carlo simulations and influence of Long-range 
interactions 

We also present some preliminary Monte Carlo (MC) sim- 
ulations that tend to confirm our RG results. We first note 
that one has to be very careful when doing MC simulations on 



a disordered system, in particular, one has to use algorithms 
that properly sample the configurations of the system, such as 
Replica MC or parallel tempering^"^. Unfortunately, these al- 
gorithms lose their full efficiency because the coefficients in 
our Hamiltonian are temperature-dependent. We thus use a 
simplified simulated annealing algorithm, which we expect to 
yield reasonable results at least for small disorder. The algo- 
rithm is the analog of the steepest descent method used in the 
literature to minimize disordered GL functionals. We empha- 
size here that such minimization methods are often inadequate 
when dealing with spin-glass-like systems as the free-energy 
typically possesses many metastable minima. We first de- 
scribe microstructure obtained in the presence of long-range 
interactions. Fig|3](a) shows typical microstructures obtained 
in different regions of the phase diagram on a 256 x 256 lattice, 
with the strength of the long-range term of Ai =4. These tex- 
tures are fully consistent with what is observed in continuum 
GL theories and in experiments. In particular, we find the 
usual cross-hatched pattern for the tweed phase. Neverthe- 
less, our RG approach provides a clear meaning to the tweed 
phase, even in the absence of long-range interactions. We 
also show in Fig. [3] (b) Field-cooling/Zero-field-cooling re- 
sults (FC/ZFC) that are usually used in both numerical studies 
and experiments to test for breaking of ergodicity and glassi- 
ness. These curves were obtained by averaging over 10^ dis- 
ordered configurations on 64 x 64 lattices, other lattice sizes 
were tested without much difference. The curves shown in 
Fig.[3](b) may be interpreted as a signature of history depen- 
dence or ergodicity breaking. However, as argued previously, 
we do not expect the 2D version of the Hamiltonian ^ to have 
a spin-glass phase. The word "spin-glass" should be under- 
stood here in the technical sense of the term, this does not pre- 
vent the system from showing kinematic "glassy behaviors" 
in the sense sometimes used by experimentalists. Therefore, 
our ZFC/FC results should be interpreted as a pure kinematic 
effect, possibly related to the slow convergence of our MC al- 
gorithm. Similar problems should occur when minimizing GL 
disordered functionals thanks to naive steepest descent algo- 
rithms. It is worth emphasizing that probing spin-glass phases 
in usual spin models is a very subtle question to address nu- 
merically, and this cannot be answered using simple ZFC/FC 
experiments. 

Finally, we suggest that the main features of our phase 
diagram persist even in the presenc e of long-range interac- 
tions (see related discussions in RefJ^^I^. The influence of 
the elastic long-range interaction on the austenite/tweed and 
tweed/strain glass transition temperatures is shown schemat- 
ically in Fig. [3];3). For no disorder (crj = 0), the austen- 
ite/martensite transition temperature decreases linearly with 
Ai, as included phenomenologically within Landau theory. 
All the transition temperatures decrease with Ai, in particular 
the glass transition is shifted to lower temperatures because 
the long range interactions compete with the randomnesJ^. 
In the asymptotic limit Ai oo, the disorder becomes ir- 
relevant and only the austenite phase remains; we therefore 
conjecture that the phase diagram is shifted to lower tempera- 
tures with increasing Ai. This result could have implications 
in the study of colossal magneto-resistance (CMR) materials 
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FIG. 2. Comparison between theoretical and experimental phase diagram, (a) Phase diagram in the temperature-disorder {r,aj) plane for our 
spin model, obtained within the RG projection approximation, r is the normalized temperature and aj characterizes the amount of quenched 
disorder in the system, (b) Experimental phase diagram of the ternary ferroelastic TisoCPdso-^Cra;). 



where the interplay of disorder and long-range strain mediated 
interactions has a bearing on phase separation of coexisting 
insulating and conducting phased. 



F. Extension to three dimensions 

We now put our work focused on the 2D SR transition in 
a more general context and explain how one can extend the 
results to more realistic situations. Although our RG analy- 
sis is in 2D and is for a square to rectangle transition (SR), 
we do not expect the salient features of the phase diagram to 
change for other transitions in 2D or 3D. As an example, let 
us discuss how one can extend our approach to the 3D Cubic 
to Tetragonal transition (CT). A spin model for this transi- 
tion was already proposed in Ref.^^. The CT transition is de- 
scribed in terms of a two-dimensional order parameter (OP) 
given by the adjunction of both deviatoric and shear distor- 



tions (e2 = (e^ 



The free energy functional can be written as a sum of 3 terms, 

F = /dV[/^(e2,e3) + /G(Ve2,Ve3)- 



fLR{e2,es] 

2n 



The Ginzburg term reads fc = ^ ( | Ve2 1 + | Ve3 1 ) , whereas 
/Li?(e2, 63) is a long-range part that shall not be important 
here (its explicit expression can be found in Ref.^^). The Lan- 
dau part is slightly more complicated /L(e2, 63) = r(e2^ + 
63 ^ ) - 2 (63 ^ - 363 62 ^ ) + (e2 ^ + 63 ^ ) ^ . This free energy can be 
minimized with respect to 62 and 63, leading to four minima 
for r < |. Let alone the 62 = 63 = minimum, the three 
other ones are given in the complex plane by e{r)ujs^, with 
1, and where similarly to the 2D case, one introduces 



s{r) 




(9) 



One can observe that r = 4/3 now corresponds to the upper 
spinodal. Retaining the Landau minima in the free energy /l, 



we define a pseudospin S such that e 
with 



Se 




(10) 



The Landau part of the free energy hence reduces to = 

(r ) QL {r)S^ (r ) , where (r ) = r - 1 + (e^ (r ) - l) ^ This 
leads to the pseudospin model 

I3H = -J{t) ^ Sj^ + A(r) ^ 5f + ^ ^ U^^Sj^, 

(11) 

With A(t) = ^(5l(t) + ee), J{r) = Do{T)e, 
and Do{r) defined as before. This model is a three- 
dimensional clock model with long-range interactions. The 
two-dimensional spin S can take 3 values on the unit circle in 
addition of the value S = 0. 

It is straightforward to extend the renormalization group 
method to this class of models. This is done by consider- 
ing hierarchical lattices with fractal dimension close to 3, for 
which the RG decimation step becomes exact if one intro- 
duces an additional coupling —K"^^. .s^{Si.Sj)'^ . We solved 
exactly the pure model defined by eq. (iT) for Ai = 0, on a 
three-dimensional hierarchical lattice to find two phases sep- 
arated by a first order phase transition around r 1, as ex- 
pected. The high-temperature austenite phase is characterized 
hy (S) = 0, whereas the martensite phase shows one of the 
three variants lying on the circle l^*! = 1. Quench disorder 
can be easily introduced through random Jij couplings and 
the RG iteration can be generalized to the disordered case. 
The calculated phase diagram is very similar to that obtained 
for SR in 2D in Fig.[2ja), although the numerical resolution of 
the RG equations is much more difficult in this case. We find 
two different "paraelastic" phases, the usual austenite charac- 
terized by the RG fixed point (J* = 0, A* = 00), and a tweed 
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Temperature r 



FIG. 3. Monte Carlo resuts, parameters are given in the text.(l) Typical microstructures obtained for Ai = 4 on a 256 x 256 lattice in the 
different phases of the phase diagram, (la) austenite, (lb) tweed, (Ic) martensite, (Id) "strain glass". (2) Example of FC/ZFC curves with 
Ai = for disorder aj=\.5, 2, 2.25, 2.5 from top to bottom. The curves represent the (normalized) susceptibility x = against the 
temperature r. (3) Qualitative phase diagram showing the influence of the long range interaction and disorder on the various phase transitions. 
Four different phases are shown: austenite, martensite, tweed and strain glass. 



fixed point with (J* = 0, A* = — oo) which corresponds 
to a disordered phase of martensitic variants. We also expect 
this kind of model to show a spin/strain glass phase, at least 
in high dimensions. These results are consistent with our ex- 
tensive study of the SR 2D model. It is thus very tempting to 
conjecture that the behavior of even more complicated clock 
models, describing more evolved transitions, should have a 
phase diagram very similar to that of the SR model. 



IV. EXPERIMENTAL PHASE DIAGRAM 

To check the predictions of our RG calculations, we ex- 
perimentally investigated the phase transformation behav- 
ior as function of temperature T and concentration x for 
Ti5o(Pd5o-a^Cra,) alloys. 

In the low Cr-content regime (x < 8), the system undergoes 
a B2^B19 martensitic transformation. The transformation 
properties of Ti5o(Pd45Cr5) alloy are shown in Figs. |4jal)- 
|4|a3). The B2^B19 martensitic transformation at its trans- 
formation temperature Mg is accompanied by a sharp in- 
crease in electrical resistivity [Fig.|4jal)]. It is found that the 
electrical-resistivity curve deviates from linearity above Ms, 
and the onset temperature of the deviation Tnd is defined as 
the tweed formation temperature. The martensitic transfor- 
mation is also characterized by a frequency-independent peak 
in internal friction [Fig. |4ja2)] and a frequency-independent 
dip in storage modulus [Fig.[4ja3)]. Note that the frequency 
independence in the Dynamical Mechanical Analysis (DMA) 
dip/peak is an important feature of martensitic transformation. 

For high Cr-content (x > 12), the alloys adopt a transi- 
tion path from austenite through tweed to strain glass. The 
transformation properties of Ti5o(Pd38Cri2) alloy are shown 



in Figs. |4jcl)- shown in Fig. |4jcl), the electri- 
cal resistivity also shows a deviation from linearity below 
Tnd , indicating the appearance of tweed state. The al- 
loy undergoes a frequency dependent storage modulus dip 
and an internal friction peak [Figs. |4jc2) and (c3)], in con- 
trast with the frequency independent behaviors during the 
martensitic transformation [Fig. |4jal)]. This demonstrates 
a dynamic freezing transition strain glass transition occurs 
in high Cr-content alloys. The ideal glass frozen tempera- 
ture (To) was obtained from by fitting the DMA dip temper- 
ature Tg{uj) and frequency uj with the Vogel-Fulcher relation 
uj = ujQ exp [—Ea/kB{Tg{uj) — Tq)]. Note that the frequency 
dependence in the DMA dip/peak is an important feature of 
strain glass transformation. 

For the crossover regime (9 < x < 12) between marten- 
site and strain glass, the alloys experience all four strain states 
upon cooling and the spontaneous transformation from strain 
glass to martensite phase (9R) takes place. Figs.|4jbl)-4(b3) 
show the predicted spontaneous transformation behavior for 
Ti5o(Pd4oCrio) alloy. Upon cooling, the tweed first form at 
Tnd = 531i^ [Fig. |4jbl)]. Further cooling gives rise to a 
frequency dispersive internal friction peak and a storage mod- 
ulus dip in the DMA results in Figs. |4jb2) and (b3), which 
correspond to a strain glass transition with frozen tempera- 
ture To of 2bOK. With a further decrease in temperature, a 
frequency-independent internal friction peak appears in the 
DMA results [Fig.[4]^b2)], which shows the similar feature as 
the martensitic transformation [Fig. |4ja2)]. This indicates a 
certain phase transformation occurs. We further studied the 
thermal hysteresis of the two phase transformations to charac- 
terize them. As shown in Fig. [5] the glass transition at higher 
temperature is associated with nearly zero thermal hysteresis 
in the internal friction anomaly and storage modulus, which 
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FIG. 4. The transformation behavior of low Cr-content Ti5o(Pd45Cr5) alloy (al)(a2)(a3), intermediate Cr-content Ti5o(Pd4oCrio) alloy 
(bl)(b2)(b3) and high Cr-content TisoCPdssCrio) alloy (cl)(c2)(c3), by means of electrical resistivity, DMA. 
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FIG. 5. DMA results on step cooling and step heating processes of 
Ti5o(Pd4oCrio) alloy, under frequency of IHz. 



is consistent with previous studies on the strain glass tran- 
sitioiPEH, From the peak interval of internal friction and 
storage modulus in Fig. [5] we found the lower temperature 
transition are associated with a thermal hysteresis about lOK, 
which shares the same feature with a first-order transition. We 
interpret the transition as a spontaneous transformation from 




FIG. 6. In-situ sychrotron XRD patterns of (a) Ti5o(Pd4oCrio) alloy 
and (b) TisoCPdssCns) alloy from 400K to 125i^^. 



strain glass to normal martensite transformation. We also em- 
ployed in-situ synchrotron XRD to check the predicted spon- 
taneous behavior for Ti5o(Pd4oCrio), as shown in Fig. |6ja). 
In the tweed region (400K and 300K), only a sharp (110) b2 
peak appears without peak splitting. At 2b2K where strain 
glass is frozen, clear peak splitting can be observed and the 
new peaks can be indexed as (002)9^^,(200)9^^ and (111)97^, 
suggesting presence of short range order associated with the 
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FIG. 7. In-situ TEM observation of Ti5o(Pd4oCrio) alloy from 330K 
to llOK. 



9R structure. With further decrease to the peak split- 

ting becomes stronger and the 9R peak height drastically in- 
creases, indicating the formation of 9R martensite. To exclude 
the possibility of the above peak splitting is caused by nan- 
odomain growth, contrast in-situ synchrotron XRD was also 
carried our on Ti5o(Pd38Cri2), in which the strain glass state 
is stable, but the system can also show the nanodomain growth 
with the decrease of temperature^^ Fig.[6];b) clearly shows 
that the (110) b2 peak keeps no splitting throughout the mea- 
surement temperature range. Therefore, the peak splitting in 
Ti5o(Pd4oCrio) alloy should indicate a transformation from 
strain glass to the martensite. 

The microstructure evolution in Ti5o(Pd4oCrio)was fur- 
ther investigated by in-situ TEM observations from 330i^ to 
llOK, shown in ¥ig.^ In the tweed region(330ir and 300K), 
contrast from nano-domains can be seen with diffused su- 
per spots in the diffraction pattern; however, when cooling to 
llOi^T (martensite region), at which 9R peaks appears in the 
Bragg reflections, parallel martensite domains are visible in 
addition to nanodomains. Erom the present experimental re- 
sults, we conclude that the second-step transition is the spon- 
taneous transition from strain glass to normal martensite phase 
9R. We note that although the spontaneous transition does oc- 
cur, the parent peak still exists in the synchrotron XRD pat- 
tern. This suggests that not all the nano-domains in strain 
glass spontaneously transform into 9R martensite within our 
measurement window. The lack of completeness of such a 
spontaneous transition can be ascribed to kinetic limitations. 
In addition, it is also noted that the spontaneous phase tran- 
sition here is rather weak and sluggish as the phase mixture 
of strain glass and martensite exists over a wide temperature 
range. This could be the reason for the absence of anomaly in 
our conventional DSC measurement (in which only latent heat 
is measured). Changes in heat capacity may aid in identifying 
the spontaneous transition. 

According to the above experimental results, a modification 
is made to the previous phase diagram of Ti5o(Pd5o-icCra^) al- 
loy s^^, where a crossover composition regime is included, as 
shown in Fig. [2] Similar phase diagram including spontaneous 



phase transition, can be found in Ti5o(Ni5o-a;Feaj) strain glass 
alloy and also in La-modified Pb(Zro. 65 Tio. 35)03 ferroelec- 
tric relaxor ceramicP^. 



V. CONCLUSION 

Our study emphasizes the importance of statistical mechan- 
ics and spin glass theory which, in conjunction with experi- 
ments, provide a general framework to understand universal 
features of the strain glass and tweed phases in ferroelastic 
materials. Our specific aim has been to emphasize how pseudo 
spin models of martensites provide a predictive route towards 
understanding aspects of glass behavior seen in experiments. 
This suggests that ferroelastics are very similar to ferroelec- 
tric and ferromagnetic materials, in the sense that they can be 
described within the same framework of statistical mechanics 
of spin models. A crucial feature of ferroelastic spin models 
is the additional S = state, which allows for two different 
"paraelastic" phases, austenite and tweed. Tweed is charac- 
terized in our study as a disordered phase of martensitic vari- 
ants. This tends to show that the tweed state observed is a 
true thermodynamical phase, that exists even in the absence 
of long-range interactions, just like martensite and austenite. 
This is an important outcome of our study as so far tweed 
has been mostly interpreted as a spin glass or as a static pre- 
cursor. We also believe that many key ideas emanating from 
the spin glass/statistical mechanics community should apply 
to the case of ferroelastics as well. For example, it still seems 
widely believed in the strain glass community that the precise 
form of the quenched disorder, modeling impurities, is crucial 
in understanding the physics of the systems under scrutiny. 
We know from experience that this is not true for spin glasses, 
and it seems highly unlikely that it could be true for disordered 
ferroelastics. We also wish to emphasize that long-range in- 
teractions, though important for microstructures, do not seem 
to be relevant to understand the global topology of the phase 
diagram. 

Although we have mainly focused on a specific model in 



2D, we argued in Sec. |III F| that our conclusions should ap- 
ply to a wide variety of 2D and 3D ferroelastic transitions; 
in particular, we believe that our calculated phase diagram is 
"generic". To be more precise, we expect to find very simi- 
lar topology of phase diagrams for other transitions with more 
variants in both two and three spatial dimensions. The alloy 
we chose as an example in the paper gives rise to two product 
phases as a function of disorder (B2 to B19, and B2 to 9R). 
Our analysis predicts that ferroelastics undergoing transitions 
to one product phase, such as TiNiFe, FePd or CaTi03, will 
show a very similar phase diagram and a spontaneous transi- 
tion. We thus expect the general topology of the phase dia- 
gram shown in Fig.jSj^S) to be quite robust. 
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